CROSS REFERENCE TO RELATED APPLICATIONS 

This application claims the benefit of U.S. Provisional Application Serial No. 60/166,242 filed 
on November 18, 1999, the contents of which is incorporated by reference as if fully set forth 
herein. 

5 BACKGROUND OF THE INVENTION 

The present invention pertains to fabrication of integrated circuits and more particularly to 
methods and systems for analyzing the performance and manufacturability of integrated circuits. 

]3 Miniaturization and scaling of integrated circuits have resulted in two conflicting trends. On one 

l; hand, margins available in high-performance designs are shrinking, increasing the likelihood that 

1 0 m inevitable variations during manufacturing will cause performance violations. On the other hand, 

s both the amount of variation as a fraction of the feature sizes, and the sensitivity of the 

^3 characteristics of transistors and interconnects to manufacturing variations is increasing. This 

li makes it essential to model and simulate accurately and efficiently the impact of manufacturing 
variations on the performance of integrated circuits. 

15 Manufacturing variations include intra-die variations, hereinafter referred to as "local mismatch", 
and inter-die variations, hereinafter referred to as "global variations". Such intra-die and global 
variations are random and systematic variations in the material composition, or in the processing 
steps that are used to fabricate the integrated circuits. Local mismatch, or intra-die variations, are 
differences in the electrical properties of circuit elements that affect components of the circuits 

2 0 that are fabricated on the same die. For example, two identically drawn transistors fabricated 



next to each other in any particular chip will still present different threshold voltage values. 
Global variations, or inter-die variations, affect integrated circuit devices by causing random 
differences in the electrical properties of circuit elements, such as transistors, resistors, 
capacitors, memory cells, wires, etc., that are fabricated on different chips from the same wafer, 
5 on different wafers, or on different batches of wafers. As circuits are designed assuming i) a 
given nominal value of the electrical properties of its elements and ii) that identically drawn 
elements behave exactly in the same way, both global variations and local mismatch affect the 
actual circuit performance and are essential to determine the yield of the product. 

H Local mismatch has not been adequately modeled in the past. Accurate matching of the electrical 

S.j 

10 m properties of active and passive elements is fundamental for functional and parametric 

f Lj performances of analog and mixed-signal integrated performance (IP) blocks, such as operational 

= & amplifiers (OPAMPs), digital to analog (D/A) and analog to digital (A/D) converters, phase 

It locked loops (PLLs), etc. Achieving target functional and parametric yield of analog and mixed- 

}t signal components frequently represents a major bottleneck for the global time-to-volume 

15 performance of complex very large scale integrated (VLSI) systems. 

Previous work in modeling and analysis of matching properties of electronic devices, such as that 
described in the paper entitled "Matching Properties ofMOS Transistors ", Pelgrom M., 
Dunjnmayer, A, and WelbersA, IEEE Journal of Solid State Circuits, Vol SC-24, pp 1433- 
1440, October 1989, which paper is incorporated by reference in this detailed description as if 
2 0 fully set forth herein, aims at deriving a suitable model of intra-die metal oxide semiconductor 
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field effect transistor (MOSFET) variance as a function of device size, layout distance and 
orientation. Although these models have become popular, their direct application in the context 
of electrical circuit simulation is difficult for two main reasons. 

First, they model the matching properties of MOSFET "macro" characteristics, such as threshold 
voltage (V Th) or saturation current (I DSAT ) that are only indirectly related to the actual "low-level" 
parameters of most widely used compact SPICE (System Program for Integrated Circuits 
Emphasis) simulation models, such as BSIM3v3 ("MOSFET Modeling andBSIM User Guide, 
Cheng, Y. andHu, C, Kluwer Academic Publishers, Boston, 1999, incorporated by reference as 
if fully set forth herein) or MOS9 ( "Compact Modeling for Analogue Circuit Simulation", 
Velghe, R., et al, IEDM Tech. Digest, pp. 485-488, 1993, incorporated by reference as if fully set 
forth herein). Therefore, a non-trivial inverse modeling process must be applied to extract the 
proper covariance structure of low-level SPICE model parameters corresponding to the available 
matching characterization data for these macro parameters. 

Second, applying a device level mismatch model to the statistical simulation of electronic circuits 
15 requires the assumption that every matched device is described by a different set of low-level 
device parameters, each associated with a corresponding random variable (RV) ("Statistical 
Modeling of Device Mismatch for Analog Integrated Circuits", Michael, C. and Ismail, M., IEEE 
Journal of Solid-State Circuits, Vol 27, No 2, February 1992; "Applying a Submicron Mismatch 
Model to Practical IC Design", Guardiani, C et al, IEEE CICC Conference, May 1994; 
2 0 "Hierarchical Statistical Circuit Characterization of Mixed-Signal Circuits using Behavioral 
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Modeling", Felt et al, IEEE- ACM International Conference on Computer Aided Design, San 
Jose, CA, November 1996, all of which are incorporated by reference as if fully set forth herein). 
The variance of the relevant circuit performance parameters can then be estimated via Monte- 
Carlo analysis. This process requires the generation of a sequence of correlated vectors of 
5 random numbers, and the evaluation of the circuit performance corresponding to each random 

vector instance either by directly using SPICE or via RSM macro-modeling. The dimensionality 
of the corresponding RV space can be very large when the simultaneous variation of all matched 
n-tuples of devices is considered. 



m The problem associated with the large dimensionality of the mismatch simulation task has not yet 

lOlrl been properly addressed. The O-space approach of Michael and Ismail (cited above), which can 

1 4 be proven to be equivalent to the Choleski factorization technique used by Felt et al. (cited 

L above), requires j£ (Nfntj) - 1) X N(pj) different RVs, where, N(m } ) is the number of matched 

ry devices of type;, N(pj) is the number of independent process factors used in the model of they th 

C3 device type, and nd is the number of different devices in the circuit. 



15 The empirical approach of Guardiani et al. (cited above), has an even greater complexity, and can 
be only applied to very simple circuits. Conti {"Parametric Yield Formulation ofMOS ICs 
affected by mismatch effect", Conti, M, IEEE Transactions on Computer Aided Design, vol 18, 
pp. 582 -596, May 1999, incorporated by reference as if fully set forth herein) proposed a method 
based on the experimental characterization of a parametrized auto-correlation function for the 

20 relevant process parameters described as spatial stochastic processes. The autocorrelation 
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function is then used to derive a symbolic formula for the system covariance matrix as a function 
of the layout parameters. Therefore, the complexity of this methodology is also proportional to 
the same number of variables as the G-space approach, however this technique is compatible 
with the statistical simulation methodology described herein, and can be used to replace the 
5 approach of deriving the component correlation matrix using a mismatch model. 

SUMMARY OF THE INVENTION 

A method of modeling and simulating the impact of global and local variation on the 
?*l performance of integrated circuits includes the steps of estimating a representation of component 
Cfi mismatch from device performance measurements in a form suitable for circuit simulation; 

1 0 W reducing the complexity of statistical simulation by performing a first level principal component 
{ ~ or principal factor decomposition of global variation, including screening; further reducing the 
!U complexity of statistical simulation by performing a second level principal decomposition 
f y including screening for each factor retained in the step of reducing the complexity; and 
O performing statistical simulation with the joint representation of global and local variation 

1 5 obtained in the step of further reducing the complexity. 

A method of modeling and simulating the impact of global and local variation on the 
performance of integrated circuits includes the following steps. First, convert, if necessary, the 
device mismatch model at the device performance level into a model suitable for circuit 
simulation, that is, at the SPICE parameter level. Second, if the global variation is provided, 
2 0 perform a first level principal-component or principal-factor decomposition and screening to 



represent the global variation. If no global variation is provided, perform this decomposition for 
the local variation alone. Third, perform a second level principal-component or principal-factor 
decomposition and screening for each independent factor identified in the second step to 
represent the local variation and reduce its dimensionality. If no global variation is explained in 
5 the initial model before local variation is applied, this step does not apply. Fourth, perform 

statistical circuit simulation and analysis with the combined set of independent factors that result 
from the second and/or third steps to estimate the impact of global and local variation jointly 
and/or separately on the circuit of interest. 

.~ (=~ 

415 

11 BRIEF DESCRIPTION OF THE DRAWINGS 

10 n Figure 1 is a schematic of a low-voltage opamp used for illustrating the method of the present 
fy invention. 

'■3 er 
f ^ 

Figure 2 is a histogram showing the impact of mismatch on the input offset voltage of the opamp 
shown in Figure 1. 

Figure 3 is a block diagram of a digital to analog converter used for illustrating the method of the 
15 present invention on circuits involving matching of more that two components at a time. 
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Figure 4 is a graph showing the efficiency of principal component analysis followed by screening 
in reducing the complexity of statistical simulation including global and local effects. 

DETAILED DESCRIPTION OF THE INVENTION 



In accordance with the present invention, an efficient representation of both inter-die and intra- 
die variation is constructed to analyze the joint impact of these sources of process variation on a 
design. The impact of manufacturing variations on a component is typically represented by 
estimating the distribution of SPICE model parameters for that component. A statistical SPICE 
\3 model that represents both inter-die and intra-die variations has to account for two types of 
'2 correlation; that is, correlation between model parameters and correlation between matched 
10 components. The correlation between model parameters arises because most commonly used 
s SPICE models utilize non-independent parameters. The correlation between different 
%J components on the same die arises because of intra-die process variation. 



15 



Both the inter-die and intra-die correlation can be represented in a single correlation matrix. 
Given n matched components (C„...,C n ) and a SPICE model for each component with m parame- 
ters (P y ,...PJ, a straight-forward implementation would construct anrnxrim correlation matrix: 



1 Pll.12 



1 



(1) 



Pl2,ll 
... ••• ... ... 

P«m,ll — PnmMm-l) * 

where the rows and columns are Tp p representing parameter j for transistor L The element p ijkl 
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represents the correlation between Tp } and T^> v 

This representation is computationally very expensive since the matrix grows as 0(nm x nm). 
Two simplifications are realized by using principal component (PC) decomposition. The first 
simplification is the use of PCS to represent the correlation between the model parameters and 
5 the second simplification is the use of PCS to represent the correlation between the matched 
devices. 




10?i 



Inter-die Principal Component Decomposition 

Principal component representation expresses a set of correlated random variables in terms of a 
set of independent random variables. This is obtained by applying a congruence transformation 
of the form: TSI 11 = A to the covariance (or correlation) matrix S, such that T is an orthogonal 
matrix and A is diagonal matrix, its elements X u are the eigenvalues of the covariance matrix. 
An important property of principal component decomposition is that: 



2 



k /-I 



(2) 



represents the fraction of the variance explained by the first q PCs. By replacing the model 
parameters (pj^p n ) by PCS (/,,.. .,/*) in the Eq. (1) one gets ankxnk size matrix, where typically 
15 k « m. The reduction in the matrix size stems from the application of a suitable threshold 

filtering algorithm based on Eq. (2). However, further simplifications are possible because using 
PCS instead of model parameters makes the correlation matrix very sparse. Because the PCS are 
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independent, the correlation between different factors for the same component is zero, i.e., p ilm = 
0 for / * m. Moreover, since the PCS are considered to represent independent sources of 
variations, correlation between different parameters on different transistors is also taken to be 
zero., i.e., p ijld = 0 for i *k and j * L This leaves only one set of non-zero entries in the 
5 correlation matrix between the same factor for different matched components, i.e., p //j7 . 

Intra-die Principal Component Decomposition 

The correlation matrix is further simplified by using a second level of PC decomposition, replac- 
ing a PC of the original model, e.g.,^, by a linear combination of a set of second-level PCS 
;i representing the intra-die variation. The result of the second PC decomposition is a unified 
10 in representation that captures both the correlation between the model parameters and the matched 
jf y components. Moreover, the transformation maintains the correlation between the model 
3 _ parameters because the RV for each original PC is replaced by an equivalent equation with the 
H same mean and variance. The second level PC transformation adds, at worst, n new PCS, 
1% resulting in a total of n x k PCS. Since the statistical SPICE model is independent of the 
15 application, the number of PCS k is fixed, and usually k « m, where m is the number of model 
parameters. This results in a mismatch simulation method that adds only O(n) new RVs for a 
circuit with n matched components. 

An advantage of the approach of the present invention is that in practice the number of factors 
required is scarcely the n x k factors required in the worst case. This happens for two reasons. 
2 0 First, not all k PCS required to capture the correlation of the model parameters are necessary for 
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representing mismatch. Usually, a much smaller number of variables are required. Second 



because the amount of mismatch is typically very small, the correlation matrix is, in general, 



characterized by a small number of dominant eigenvalues. Therefore, using the equation that 



gives the fraction of total variance explained by the k-th PCS, the variance of the system can be 



5 



well approximated by using a small number of PCS compared to the worst-case n. 



For example, in the asymptotic case of perfectly matched devices, this method automatically pro- 
duces only one PC (corresponding to the dominant eigenvalue) for each independent process fac- 
« tor, resulting in no increase in the number of RVs. In the intermediate situations between the 
m worst-case and the asymptotic case, the approach of the present invention provides an 
10 in approximation scheme where a small number of PCS can be selected to accurately approximate 
■ y the correlation between the components. 



[3 Mismatch characterization typically does not produce the correlation between PCS of a statistical 
SPICE model Usually, the result of mismatch characterization is a set of coefficients for a mis- 
15 match model of device performance. For example, MOS transistor mismatch characterization 
often results in the coefficients of a Pelgrom-style model for threshold voltage (V^) and satura- 
tion current (l DSAT ) or the transistor gain-factor (/:') (see the Pelgrom, et al reference cited above). 
To utilize this information in the representation described above it must be converted into corre- 
lation between PCS of the statistical SPICE model. 



Statistical SPICE Models with Mismatch 
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The procedure for determining the unified statistical representation starts with a statistical SPICE 
model for correlation between model parameters. A subset of the PCS of this model is selected 
for representing component mismatch. The selection is based on two considerations: the device 
characteristics for which mismatch has been characterized, and the weightings of the different 
PCS in the equation for each SPICE parameter that impact these device characteristics. For 
example, if mismatch characterization has been performed for threshold-voltage and gain-factor 
of longchannel MOSFETs, then PCS that have the most impact on the VTHO and UO (mobility) 
parameters of the BSIM3v3 model are selected. Once the PCS of the model without mismatch 
have been selected, a numerical optimization is performed to find the correlation between these 
factors in order to obtain the measured mismatch in device characteristics. 

Each optimization results in the correlation of the selected PC for one set of device geometry and 
layout parameters. Repeating this procedure for different geometries and layout distances results 
in a set of correlation values that are either fit with an interpolating function or represented in a 
look-up table. This procedure results in a model of PC correlation as a function of device 
geometry and layout. The step of modeling PC correlation as a function of device geometry and 
layout is performed only once for a particular technology and layout style. 

The method of the percent invention can be used as part of the Circuit Surfer statistical design 
and verification environment. This environment supports many tasks typically required for 
analog and mixed-signal design for manufacturability, such as statistical simulation, sensitivity 
analysis, response surface modeling, and circuit optimization for manufacturability. Mismatch 




simulation is implemented in this environment as an annotation of the circuit netlist to specify 
the matched components, and annotations to the statistical SPICE models to include the effect of 
mismatch. The annotated netlist and SPICE models are used to derive a separate statistical 
SPICE model for each component using the two-level PCA. The modified netlist forms the input 
5 to Circuit Surfer. This implementation enables all the capabilities of statistical design for 

mismatch analyses such as variable screening, response surface modeling and Monte-Carlo using 
mismatch factors, and optimization of a design to reduce its mismatch sensitivity. 

The use of the present invention is illustrated with two applications: a low voltage OPAMP and a 
%j D/A Converter. These topologies were implemented based on their continued significance in 
10 " j low voltage VLSI signal processing applications and their performances sensitivities to device 
" mismatch characteristics. 

J5 Low Voltage OPAMP 

u The following is a description of an embodiment of the method, in accordance with the present 
u invention, for statistical modeling and simulation of the impact of global variations and local 
15 mismatch on the performance of operational amplifiers (OPAMPs) of the type described. Such 
local mismatch and global variations are typically caused by random and systematic variations in 
the material composition, or in the processing steps that are used to fabricate the integrated 
circuits. Global variations affect integrated circuit devices by causing random differences in the 
electrical properties of the circuit elements, such as transistors, resistors, capacitors, memory 
2 0 cells, wires, etc., that are fabricated on different chips from the same wafer, on different wafers, 
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or on different batches of wafers. The local mismatch, or intra-die variations, are differences in 
the electrical properties of circuit elements that affect components of circuits that are fabricated 
on the same die. For example, two identically drawn transistors, fabricated next to each other in 
any particular chip, will still present different threshold voltage values. 

5 The OPAMP is of the 2-stage rail-to-rail class-AB architecture shown in Figure 1 which employs 
a constant-g m input stage with tail current control. The commonality of complimentary input 
stages in contemporary designs makes this topology ideal for exploring the effects of mismatch 
on offset voltage (Vos) given both NMOS and PMOS mismatch effects. In addition, low voltage 
JS design examples such as this help exploit the increasing dependence of proper signal resolution 
1 0 in on low Vos. A Pelgrom model for the mismatch in the threshold voltage and gain factor of long- 
fSJ channel transistors is available for this class of technologies. For the matched NMOS pair Ml- 
^ M2, the model specifies o(A( Vth)) = 2mV and a(A(fc')) = 0.7 %. For the PMOS pair M3-M4, the 
U mismatch was: (o(A(Vta)) = LSmVand oCAC* 1 )) = 0.45 %, where Vth represents the MOSFET 
O threshold voltage and k' represents the MOSFET gain factor. 

15 A statistical SPICE model was estimated for the process technology used in this example. The 
model was estimated using measured variance and correlation between key transistor 
performances. This set of performances includes: transistor drive current in the linear and 
saturation regions (Idsat, Idlin) and transistor threshold voltage in the linear and the saturation 
region (Vtsat, Vtlin). The statistical SPICE model is represented in terms of principal factors, 

2 0 thus achieves the effect of first-level principal component decomposition. 
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The mismatch model described above is not suitable for direct use in circuit simulation. A 
SPICE level representation of mismatch is derived by selecting SPICE model parameters that 
effect long-channel threshold voltage and the gain factor. For BSIM3v3 SPICE models, these 
parameters are VTHO, TOX, and UO. Numerical optimization is performed to find the 
5 correlation between these factors in order to obtain the measured mismatch in threshold voltage 
and gain factor. 



The statistical model including mismatch was used to simulate the impact of mismatch on key 

0 performance measures on this design. This simulation makes use of the method of this invention 
J; 1 : to perform the second level principle component decomposition. The impact of mismatch on this 

10 design is shown in Table 1. 

? ; h 

1 is" 

C3 Table 1: Impact of Mismatch on OPAMP (10,000 sample Monte-Carlo) 



Performance 


Without Mismatch 


With Mismatch 


Mean 


std. dev 


Mean 


std. dev 


Vos (mV) 


-1.803 


0.09627 


-1.802 


0.7216 



A unified model is derived by applying the method described in the previous sections. The first 
15 level PCA results in nine first level PCS (pcal...pca7, pcaLint and pcaWint) for the four 

matched transistors, thus resulting in a 36 x 36 sparsified covariance matrix. A second level 
PCA is performed using this matrix. Three dominant eigenvalues explain over 98% of the total 
variance. In this way, it is possible to account for both inter-die and intra-die variations with the 
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addition of three variables. These three variables are suitably obtained by re-defining pca2, pca4 
and pca5 as linear combinations of two independent components. 

Two sets of Monte-Carlo simulations were performed, one using only the inter-die model and the 
second using the unified model. As expected, this design shows that Vos is extremely sensitive to 
5 mismatch. By quantifying the exact impact of mismatch the method of the present invention 
allows an accurate assessment of manufacturability of this topology. It also shows that 
parametric yield estimates can be overly optimistic for designs sensitive to mismatch if mismatch 
C3 effects are not accounted for in statistical simulation or worst-case models. Figure 2 shows the 
J r [ distribution of Vos with and without mismatch. The overlay illustrates the 7X increase in the 
1 0 offset voltage standard deviation. Without the ability of efficient and accurate mismatch analysis 
O the impact on parametric yield due to this increase would be missed. 

E 

SL s j: 
t Jj 

s ir 
r. 1 

zi D/A Converter 

The D/A Converter shown in Figure 3 uses a binary weighted current implementation with eight- 
bit resolution. The output current is summed via current switching and output to a linear 
15 current/voltage (I-V) converter. Data converter non-linearity is a measure of the error induced by 
the converter and it is sensitive to the matching between the transistors comprising the ratioed 
current source. 
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Simulation with and without intra-die variation was carried out as described for the OP AMP. 
The procedure was deriving the statistical SPICE model using principal-factors and estimating a 
representation of device mismatch in form of a SPICE model are identical to those for the 
OP AMP. The simulation was carried out using the method of this invention. The mismatch 
5 effects were examined for the binary weighted current source identified in Fig. 3. Specifically, 

the analog output was tested for differential non-linearity (DNL), defined as the deviation of each 
set of adjacent steps at the analog output from their ideal value (1 LSB). Here, we have defined 
the performance specification to be +/- Vi LSB. The results of statistical simulation with and 
p without intra-die effects are shown in Table 2. 

1 0 l 2 Table 2: Impact of mismatch on D/A Converter (10,000 sample Monte-Carlo) 



Performance 


Without Mismatch 


With Mismatch 


Mean 


std. dev 


Mean 


std. dev 


DNL (LSB) 


-0.176 


0.086 


-0.198 


1.155 


Yield 


100% 


34% 



Intra-die variation causes the standard deviation of DNL to increase by a factor of 13. Such a 
15 large increase in the increment causes large changes in the linearity of the D/A Converter. As in 
the previous example, including intra-die effects is vital to get accurate estimation of circuit 
performance variability, but since the circuit's performance for this topology depends upon the 
pairwise matching of eight different transistors, a complete statistical simulation using a standard 
methodology would entail running a Monte-Carlo experiment with large number of correlated 
2 0 RVs. 
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For example, a circuit with eight matched transistors and 7 to 12 independent sources of process 
variation (typical of advanced CMOS processes) requires the generation of a sequence of corre- 
lated random numbers for as many as 96 different RVs. Several thousand SPICE evaluations 
would be required to stabilize the Monte-Carlo results for such a high-dimensional system. The 
advantages of complexity reduction using second order principal component decomposition are 
illustrated by the graph depicted in Figure 4. This figure shows the standard deviation of DNL 
versus the number of additional factors employed for mismatch simulation. In fact, although in 
the worst case the present technique would require 18 additional RVs, by applying the screening 
methodology previously described under the heading "Intra-die Principal Component Decompo- 
sition" above, it is possible to model very accurately the impact of intra-die variability on the 
D/A Converter with only two additional variables. This result makes the simulation of mismatch 
applicable to an inherently larger set of matched transistors; a requirement as supply levels and 
feature sizes continually decrease and previously negligible device effects become increasingly 
significant. 
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